Absence of electron dephasing at zero temperature 
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Dephasing of electrons due to the electron-electron interaction has recently been the subject of 
a controversial debate, with different calculations yielding mutually incompatible results. In this 
paper we prove, by means of Ward identities, that in a system with time reversal invariance, neither 
a Coulomb interaction nor a short-ranged model interaction can lead to phase breaking at zero 
temperature in spatial dimensions d > 2. 
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I. INTRODUCTION 

Many interesting phenomena in electronic systems at 
low temperatures depend on the electrons' phase coher- 
ence. Among these are the 'weak-localization effects',^ 
which are nonanalytic dependences of zero-temperature 
(T = 0) electronic correlation functions on frequency 
and/or wavenumber. Loss of phase coherence cuts off 
these singularities, and transfers the nonanalyticities into 
nonanalytic dependences on the phase breaking parame- 
ter. The standard interpretation of the weak- localization 
experiments is based on the temperature dependence of 
various phase breaking processes and the resulting non- 
analytic temperature dependence of various observables,^ 
and an understanding of phase breaking processes is 
therefore of great importance. Of particular interest is 
phase breaking due to inelastic scattering caused by the 
electron-electron interaction, since at low temperatures 
this dominates other temperature dependent scattering 
processes, e.g., those caused by phonons. 

For many years there had been consensus that at zero 
temperature the phase-breaking rate due to interactions 
vanishes, so that in the absence of other dephasing pro- 
cesses like, e.g., magnetic scattering, the electrons main- 
tain full interference capabilities. This consensus was 
based on calculations of various phase-breaking times. 
An important example is the single-particle inelastic life 
time Tin , which describes the decay of the single-particle 
Green function, 
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Here p is the wavenumber, — e(p) — fi with e(p) the 
renormalized single-particle energy and fi the chemical 
potential, Un = 27rT(n -I- 1/2) is a fermionic Matsubara 
frequency, and Z is the wavefunction renormalization. 
Landau Fermi-liquid theory predicts the single-particle 



inelastic scattering rate, l/rin, in clean electron fluids to 
vanish like l/rin oc T^.^ Schmid calculated 1/rin for dis- 
ordered systems with diffusive electron dynamics, and 
found it to be enhanced compared to the clean case, 
vanishing like T'^/-^ in 2 < d < 4 spatial dimensions.^ 
The transport properties measured in weak-localization 
experiments are governed by a similar relaxation time, 
namely, the phase breaking time that is related to a 
two-particle correlation function often referred to as the 
diffuson propagator, ^'^ 
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Here fin — 27rTn > is a bosonic Matsubara frequency, 
if is a frequency renormalization, and D is a renormal- 
ized diffusion coefficient. For 2 < d < 4, many different 
calculation found the temperature dependences of and 
Tin to be qualitatively the same, 1/r^ oc l/ri„ oc T'^/'^J In 
a weak- localization context, instead of the diffuson, the 
closely related Cooperon propagator is often discussed, 
which has a mass that in time-reversal invariant systems 
is also given by t^.* An overview within a unifying the- 
oretical framework has been given in Ref. 9. 

This consensus was shattered in the wake of experi- 
ments on disordered metallic wires that found that 1 /r,^ 
saturates as T — > 0.^"'^^ The experimentalists assert to 
have ruled out other sources of dephasing, which suggests 
that the electron-electron interaction can lead to resid- 
ual dephasing even at T = 0. These experiments deal 
with finite-size one-dimensional or two-dimensional sam- 
ples only, and for these systems an explanation in terms 
of quantum fluctuations has been proposed. However, 
more detailed theoretical work along the same lines found 
such a residual contribution to l/r^ not just in c? = 1 and 
d — 2, but also in d = 3,^^ and all of these results are in 
contradiction to the earlier theoretical work.® The quan- 
tum fluctuations that are proposed to lead to the satu- 
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ration have been taken into account in the earher calcu- 
lations, as Rcf. 9 makes explicit, yet the results disagree. 
This contradiction has recently given rise to an intense 
debate. ^"^'^^ It is important to note that this disagreement 
by far transcends the technical issue of what the behav- 
ior of a particular observable is, and how to calculate 
it. A finite dephasing rate at T = in three-dimensions 
would mean the breakdown of a Fermi-liquid picture for 
disordered bulk Fermi systems, with far-reaching conse- 
quences. 

In this paper we attack the problem from a very gen- 
eral theoretical angle. Rather than calculating the phase- 
breaking rate explicitly, we use general symmetry argu- 
ments and Ward identities to prove that the electron- 
electron interaction cannot give rise to a nonvanishing 
value of l/r^ at zero temperature in d > 2. This mathe- 
matical result establishes that the diffuson propagator is 
always massless at zero temperature. We then show that 
in the presence of time reversal invariance, the same is 
true for the Cooperon propagator. 

II. SPONTANEOUS SYMMETRY BREAKING 
AND SOFT MODES IN ELECTRON SYSTEMS 

The discussion in the literature of conflicting results 
obtained by different calculations that we summarized 
in the Introduction fails to mention an important point. 
Namely, the question of the existence or otherwise of a 
finite electron dephasing time at zero temperature is a 
question about the soft-mode spectrum of the system, 
since a nonzero dephasing rate gives certain propagators 
a mass, see Eq. (1.2). The existence and identity of soft 
modes, however, is governed by fundamental properties 
of the system and can be studied rigorously, and in very 
general terms, without reference to any particular ex- 
plicit calculation. In fact, such an analysis should always 
be performed before carrying out calculations, since the 
former provides important checks for the latter. This is 
particularly important if the calculations involve approx- 
imations, as it is crucial to make sure that the approx- 
imations are 'conserving' in the sense that they do not 
violate any exact qualitative properties. Of course, sys- 
tematic perturbation theory, if done correctly, is always 
conserving in this sense, as any exact properties will hold 
order by order in perturbation theory. However, given 
the complexity of many perturbative calculations, and 
the fact that perturbation theory by itself can never es- 
tablish that some property is actually exact, it is still 
valuable to have exact benchmarks. 

There are two known mechanism for guaranteeing that 
certain correlation functions are massless, viz., (1) con- 
servations laws, and (2) spontaneously broken symme- 
tries. Conservation laws ensure that a conserved quan- 
tity cannot be locally created or destroyed. An excitation 
that couples to a conserved quantity can only decay by 
some of the quantity being transported into or out of the 



perturbed region. For long-wavelength excitations, this 
takes a long time. Therefore, the correponding correla- 
tion functions will decay infinitely slowly in the limit of 
zero wavenumber and zero frequency, i.e., they are mass- 
less. An example are density fluctuations in any isolated 
many-particle system, whose decay rate goes to zero with 
vanishing wavenumber due to the particle number con- 
servation law. This is of course true both at zero and 
at nonzero temperature. A spontaneously broken con- 
tinuous symmetry also leads to soft modes, the number 
of which depends on the number of parameters of the 
symmetry group that is broken, and on the extent of the 
breaking. This is the content of Goldstone's theorem, 
and the resulting soft modes are usually referred to as 
Goldstone modes. An example are the spin waves in 
a Heisenberg antiferromagnet, which are the Goldstone 
modes resulting from the spontaneously broken spin ro- 
tation symmetry. 

The systematic analysis of symmetries and the re- 
sulting soft modes in disordered itinerant electron sys- 
tems was pioneered by Wegner,^'' and by McKane and 
Stone. These authors showed that, for noninteracting 
electrons, the diffusive excitations that give rise to the 
weak-localization effects are the Goldstone modes of a 
spontaneously broken symmetry that can be formulated 
as a rotation in frequency space. Roughly speaking, it 
is the symmetry between retarded and advanced Green 
functions that is spontaneously broken whenever there 
is a nonzero density of states. A nonzero temperature 
serves as an external field conjugate to the order param- 
eter, i.e. the density of states, and gives the Goldstone 
modes a mass. The results of Refs. 13, 14 amount to the 
claim that either, the electron-electron interaction pro- 
vides such a conjugate field even at T = 0, making the 
Goldstone modes massive, or, that time reversal invari- 
ance is always broken in interacting many-body systems, 
making the Cooperon different from the diffuson. 

Wegner's symmetry analysis has been generalized to 
the case of interacting electrons in Refs. 19, 20, and was 
used to construct an effective field theory for disordered 
interacting electrons. Although it is implicit in this work 
that at zero temperature, and in d > 2, the electron- 
electron interaction cannot lead to any dephasing, the 
controversy discussed above makes it worthwhile to re- 
analyze the symmetry arguments and prove this state- 
ment explicitly. We will now proceed to do so. 

III. ABSENCE OF PHASE BREAKING AT ZERO 
TEMPERATURE 

A. Field theory for fermions 

Our starting point is a general field theory for elec- 
trons. For any fermionic system, the partition function 
can be written^^ 
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z = J exp(5[v;,v]) , (3.1a) 

where S is the action in terms of the fermionic (i.e., 
Grassmann valued) fields -ip and ip. We consider an ac- 
tion that consists of a frcc-fermion part Sq, a part Sdis 
describing the interaction of the electrons with quenched 
disorder, and a part Sint describing the electron-electron 
interaction. 



S — So + Sdis + Sint 



(3.1b) 



Each field ip or ip carries a Matsubara frequency index 
n and a spin index a =t, [■ In terms of these fields, 
reads explicitly 



So 



T 

(x) [iu)n + e(V) + IJ] tpnai^) 



(3.2) 



with n the chemical potential, and a;„ = 27rT(n + 1/2) 
a fermionic Matsubara frequency, e denotes the disper- 
sion relation. For instance, for free electrons one has 
e(V) = — V^/2m, with m the electron mass. Since we 
will be interested in four-fermion correlation functions, 
it is useful to introduce a matrix of bilinear products of 
the fermion fields, 



B„ 



/-iJjnlAn^ -^^nlAnl -%/^'nlAnl V^'nT^-'mt \ 

2 V^ni-^mT iinlAnl 'ipnl'^"^l ""^"-iV'mT 

\--0„TV'mT -V'ntV'mi -V'ntV'mi V'ntV'mT / 

Qnrn ■ (3-3) 



where all fields are understood to be taken at position x. 
The matrix elements of B commute with one another, 
and arc therefore isomorphic to classical or number- 
valued fields that we denote by Q.^^ This isomorphism 
maps the adjoint operation on products of fermion fields, 
which is denoted above by an overbar, on the complex 
conjugation of the classical fields. We use the isomor- 
phism to constrain B to the classical field Q, and exactly 
rewrite the partition function^^ 

Z = J £»[VS,V] e^[^'^] J D[Q]S[Q-B] 

= y£>[VS,V] e^I'^'^J J D[Q]D[A] e^[^(«-^)l 

= j D[Q] D[A] e-^[«'^] . (3.4) 

A is an auxiliary bosonic matrix field that serves to en- 
force the functional delta-constraint in the first line of 
Eq. (3.4), and the last line defines the action A. The 
matrix elements of both Q and A are spin-quaternions 
(i.e., elements of Q x Q with Q the quaternion field). 
Prom Eq. (3.3) we see that expectation values of the 



Q matrix elements yield local Green functions, and Q- 
Q correlation functions describe four-fermion correlation 
functions. The physical meaning of A is that its expec- 
tation value plays the role of a self energy, see Ref. 19. 

It is convenient to expand the 4x4 matrix in Eq. (3.3) 
in a spin-quaternion basis, 



Q 



nm 



(x) = E (r. 

r,i=0,3 



) Si) j.Qnm (^) 



(3.5) 



and analogously for A. Here tq = sq = I2 is the 2x2 unit 
matrix, and tj = —sj = {j = 1,2,3), with (Ji.2,3 

the Pauli matrices. In this basis, i = and i = 1,2,3 
describe the spin singlet and the spin triplet, respectively. 
An explicit calculation reveals that r = 0, 3 corresponds 
to the particle- hole channel (i.e., products '!/'f/')> while r = 
1,2 describes the particle-particle channel (i.e., products 
iIjiIj or We will be particularly interested in the 

matrix elements of qQ and 5*5, for which the isomorphism 
expressed in Eq. (3.3) reads 

0012 (x) - ^ X! [V'la(x)V'2a(x) + ^2a{y^)-4^1a{y^)] , 

(T 

(3.6a) 

i<3i2(x) = \ [V'iT(x)'02i(x) -I- Vi2T(x)V;u(x) 

+ ViT (x) V2i (x) + V2T (x) Vu (x)] . (3.6b) 

Prom the structure of Eq. (3.3) one obtains the follow- 
ing formal symmetry properties of the Q matrices, 



(3.7a) 
(3.7b) 
(3.7c) 
(3.7d) 
(3.7e) 



Here the star in Eq. (3.7e) denotes complex conjugation. 

By using the delta constraint in Eq. (3.4) to rewrite 
all terms that are quartic in the fermion field in terms of 
Q, we can achieve an integrand that is bilinear in ij) and 
■0. The Grassmannian integral can then be performed 
exactly, and we obtain for the action A 

A[Q, A] = ^dis + Ant + ^ In (<^o ^ - 

+ j dxtT (a(x) Q(x)) . (3.8a) 



r^nrn 
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(r = 0,3; i 


= 1,2,3) 


r^nm 


~ rQmn ' ('^ ^ 


1,2) , 




r^nm 


~ ~rQmn ' 


= 1,2; i = 


1,2,3) 


r^nm 


— -'O 

r-T! — n— 1,— m— 1 







Here 



Go' = -a, + e(V)+M 



(3.8b) 



is the inverse free electron Green operator, with dr the 
derivative with respect to imaginary time. Tr denotes a 
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trace over all degrees of freedom, including the continu- 
ous position variable, while tr is a trace over all those 
discrete indices that arc not explicitly shown. 

For explicit calculations, the electron-electron inter- 
action ^int is usually decomposed into four pieces 
that describe the interaction in the particle-hole and 
particle-particle spin-singlet and spin-triplet channels, 
respectively. For the purposes of the present paper, 
which will study exact structural properties rather than 
explicit calculations, this decomposition is neither nec- 
essary nor desirable. We therefore write the interaction 
part of the action as the basic density-density interac- 
tion between electrons mediated by a statically screened 
Coulomb potential. In terms of fermionic fields, this 
reads 

•S'int = ~^ y dxdyvsc{x-y) ^ ■^ai (x) (x) 



xV'<72(y)V'a2(y) • (3.9a) 



The statistical properties of the random hmction m(x) 
are governed by a distribution functional P[u]. The par- 
tition function, Eq. (3.4), is a functional of ■u(x), and any 
physical quantity X, like the free energy, or any corre- 
lation function, needs to be averaged over the disorder 
according to 



D[u{x)]X[u]P[u] 



(3.10c) 



Here 



In Refs. 19, 20 the disorder average was performed by 
means of the replica trick. Here, since we are not con- 
cerned about explicit calculations, we will formally per- 
form the disorder average by means of the basic defini- 
tion, Eq. (3.10c). 

The above definitions completely specify our model. 
We stress that the Q-matrices merely provide a conve- 
nient shorthand for bilinear products of fermion fields. It 
would be perfectly possible to proceed in terms of the lat- 
ter, and whether or not one introduces the matrix fields 
is a matter of taste. 



Mx) = r^e- 



-«|x| 



(3.9b) 



is the screened Coulomb potential, with e the electron 
charge and k the Thomas- Fermi screening wavenumber. 
We stress that we use a screened interaction for conve- 
nience (and since screening is a real physical effect) only, 
and that our arguments do not rely on this. A suitable 
modification of the procedure in Sees. Ill B and III C still 
applies if one works with a bare Coulomb interaction, and 
the results are the same, see Ref. 20. 

Its Fourier transform is, in d spatial dimensions, 



1 



jd-i) 



Nf [q2 -hK2](rf-l)/2 ' 



(3.9c) 



with A^F the density of states of clean, noninteracting 
electrons at the Fermi level. In terms of the Q matrices, 
Eq. (3.9a) reads 

^int = 8T / dxdy ^ J2 ^sc(x - y) 

r=0,3ni,n2,m 
XrQn.,n.+mi^)rQn.,n.+miy) > (3-9d) 

Finally, the interaction of the electrons with the static 
disorder is given by a random potential w(x) that couples 
to the electronic density, 

-^dis = - dx U{x) ^ ^ V'„,a(x) V'n,<r(x) , (3.10a) 
an 

or, in terms of the Q, 

Adis = [dx ti(x) ^ [Sg„„(x) . (3.10b) 



B. Symmetry Considerations, and a Ward identity 

A symmetry analysis of the action, Eqs. (3.8) - (3.10), 
has been given in Refs. 19, 20. Here we recapitulate 
those aspects of this analysis that are important for our 
purposes. 

Let us consider an infinitesimal transformation of the 
Q-matrices, Q TQT''^, with T = CfC'^, where C = 

iTi ® S2, and 



nni ^mn2 



(3.11) 



ril > and n2 < are fixed frequency indices that char- 
acterize the transformation. This amounts to a rotation 
in frequency space by an infinitesimal angle 9. These ro- 
tations are elements of the symplectic symmetry group 
Sp(8-/V) over the complex numbers, with 2N the num- 
ber of frequency labels, which governs the symmetry 
properties of the action. Under such a rotation, the 
Q-matrices transform like 



(3.12a) 



with 



SQnm = {^nniQn^m + SmniQnn2 — 2)) 9 + 0{9'^) . 

(3.12b) 

The A transform analogously. The symbol (1 ^ 2) de- 
notes the same terms as written previously, but with the 
indices 1 and 2 interchanged. 

Of the terms in the action, Eq. (3.8a), only ^int and 
the Tr In term are not invariant under these rotations. 
To linear order in the transformation parameter 9 one 
finds A—i' A + SA, with 
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(3.13a) 



with G = {Gq^ - iA)-\ and 

[Siuj)^l^ = {S-fini^mn2 ~^ ^nn2^mni) ^^m— "2 i (3.13b) 
JAnt = 32 I dxdy z;,e(x - y) ^ ^ E [r^n^n, W 

r=0,3 UaTli, 

x°Qn.,«.-(n„-«.)(y)-(l""2)]e . (3.13c) 



Now introduce a source J for the Q-fields, and consider 
the generating functional 



14) 



By performing an infinitesimal transformation on the Q 
and the A, one obtains from Eq. (3.14) 







= J D[Q]D[A] 



^■^ + j rfxtr(j(x) SQ{x)) 

^g^+/dxtr(j(x)Q(x)) 



(3.15a) 

Differentiating this identity with respect to oJ„4„3, > 
0, < 0, and putting J = yields 

= (m0Q„^„^(x)) + (°(5O)„^„Jx)) , (3.15b) 

Here the angular brackets denote an average with respect 
to the action A. Prom Eq. (3.13a) we see that the first 
term involves an average (GQ). By using the identity^'' 
{GQ) = —2i{Q^), and performing the disorder average, 
Eq. (3.10c), we finally obtain 

^^ni— 712 -^nin2, 713714 ^711713^712 714 "^711712 ^^7ll7l2 ,713714 

(3.16a) 



Here 



dis 



(3.16b) 



D. 



711712, 713714 



= J dy {(°(AQ)„^„Jy) 

Xo(Ag)„3„,(x))}^.^ , (3.16c) 
with AQ = Q - {Q), and 



711712 ,713714 



-32w,e(k^0)A/-„,„,T J2 ^n.- 



■nt,^ni—n2 



xD. 



-32 ^ / dydzv,,iy-z)T ^ 

r=0,3"' 7i„,7i(, 

{(°(AO)„„„,_(„„_„,)(-q):!(AQ)„^„^(q) 



xg(AQ)„3„,(x))}^^^-(l-2) 



(3.16d) 



The physical meaning of the propagator _D is a phase 
space density correlation function. The homogeneous, 
static density susceptibility x can be obtained from D 
by the relation^^ x = ^^TJ2n,m Dnn.mm- 

The Ward idcntity,i^'2° Eq. (3.16a), relates the two- 
point Q-correlation function D to the one-point function 
N and the three-point function W. Various aspects of 
this identity have been discussed in Refs. 19, 20. In the 
next subsection we focus on the implications for the de- 
phasing time. 



C. Absence of Dephasing 

We now use the Ward identity, Eq. (3.16a), to show 
that there cannot be a nonzero phase relaxation rate at 
zero temperature in > 2. We will first concentrate 
on the diffuson propagator, and then show that, in time 
reversal invariant systems, the same conclusions hold for 
the Cooperon propagator. 



1. The diffuson propagator 

Let us consider the structure of the Ward identity, Eq. 
(3.16a). Equation (3.6a) shows that (oQ„„(x)) is pro- 
portional to the on-site electron Green function, so for 
ni > 0, 712 < 0, f2„i_„2 = - 0, iV„i„2 ap- 

proaches a constant N that is proportional to the density 
of states at the Fermi level. D and W can be decomposed 
into two pieces each that have different frequency struc- 
tures, 

D = S S n(<i<:) 

-^711712,713714 — "7ll,7l3 "712,714 "'^711712 

+^711-712 ,713 -714 -^711712 ,713714 ' (^•l'^) 

and analogously for W. The superscripts '(c)' and '(dc)' 
refer to 'connected' and 'disconnected' (at a fermionic 
level) contributions. The disconnected piece of D, D'^'^'^\ 
is identical, apart from a proportionality constant, with 
the homogeneous limit (p — > 0) of the diffuson propa- 
gator L, Eq. (1.2), that was defined diagrammatically in 
Rcf. 5. For nonintcracting electrons, the connected piece 
of D vanishes, and we have, in the limit flni-n2 ~^ 0, 



(3.18) 



This expresses the undisputed fact that in the absence of 
an electron-electron interaction, the diffuson is massless 
at zero temperature.^^ 

The question is now whether the presence of W can 
change this fact. The only possibility is a cancellation 
between N and W in the limit of zero frequency and 
zero temperature. To see that this cannot happen, let us 
consider the equation for D^'^'^^ in this case. It is sufficient 
to show that _C)('^^) remains massless, since i?^'^^ and _C)('^^) 
cannot cancel one another due to their different frequency 
structures. Prom Eq. (3.16a) we obtain 
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so n(dc) = TV - W^^'^) 
An evaluation of Eq. (3.16d) yields 



(3.19a) 



W(<*'=) = 8 lim / dy dz v^^iy - z)T V 

X |g(x, Z; ) {ipna {y)so'ipn, (y)) 

X ((V;„2_„„+„,(z)soV'n2(x)) ) -(1^2) 



dis 

(3.19b) 

Here we have introduced two-component spinors, 
ipn{'x.) = {'ipn\i^),il'ni{^)), ^nd a scalar product in spinor 
space that is given by the usual matrix product. The 
{-tpip-tpipy are connected i/;-correlation functions, with the 
cumulant taken with respect to the quantum mechanical 
average only, and G(x, z;iw„) = {'tpna{^) i^nai^)) is the 
Green function for a given disorder configuration. 

Now consider W'-'^''\ Eq. (3.19b), as a function of the 
screening wave number k. The behavior for small k is de- 
termined by the infrared properties of the integral that 
defines W^'^''\ The connected four-fermion correlation 
functions contain at least one screened interaction po- 
tential, and any four-fermion function is at most diffu- 
sive, i.e. diverges at most like an inverse wavenumber 
squared for small wavenumbers. Counting the number 
of frequency and momentum integrals in Eq. (3.19b), we 
conclude that the leading zero-temperature contribution 
to T4^(<^'=) for K — > is bounded above by 



W^'^"') < const. / dpp"^-^ 



,2(d-l) 



1 



'^-^ + o{k'^-^) , (3.20) 



OC K 



where o{x) denotes terms that vanish faster than x as 
a; ^ 0. A is an ultraviolet cutoff that is on the order of 
the Fermi wavenumber. It follows that VF'^'*'^'' vanishes 
continuously as k for d > 2. The density of states, 
in contrast, approaches a finite value, and so = 0{\) 
in this limit. We see that TV and W'^'^"^ are different func- 
tions of K, so in general they cannot cancel. This means 
that the diffuson propagator remains soft in the pres- 
ence of interactions. It can become massive at most at a 
special value of the screening wavenumber, which would 
signalize a phase transition into a non-Fermi liquid state. 
If one replaces the screened Coulomb interaction by a 
different short-ranged model interaction with interaction 
amplitude F,^'' then all of the above arguments still apply 
if one considers N and H^f^'^) as functions of V instead of 

K. 

The realization that W and N arc necessarily different 
functions of the interaction, and therefore cannot gener- 
ically cancel, is the most important ingredient of our 
proof. Notice that a cancellation would require 
be of 0(1) for k ^ 0. It follows from Eq. (3.20) that the 



only way for this to happen is to have a stronger than dif- 
fusive infrared singularity in a four-fermion propagator. 
However, the nonvanishing phase relaxation rate that 
might result from this would make the diffusive Gold- 
stone modes massive, i.e., it would lead to a theory that 
is less singular in the infrared. Any scenario that has 
the diffuson become massive due to the electron-electron 
interaction would therefore not even be internally consis- 
tent. 

We finally note that no cancellation between A^ and 
^(dc) jg necessary to produce dephasing at a nonzero 
temperature. According to the relation between _D('''^) 
and L from Eq. (1.2) mentioned above, and after an an- 
alytic continuation to real frequencies f2, one has 



const. 



H + i/Qr^ 



(3.21) 



As long as, for T ^ and O — » 0, a regime exists where 
rJr^ >> 1, Eq. (3.21) allows for an expansion in powers 
of with the leading contribution to A^ — W'^'^'^^ 

a nonzero constant. This is the case if 1/r^ vanishes 
continuously as T — > 0, and in particular if 
However, in order to have — *■ const, for T ^ 0, the 
leading 0(1) contribution must cancel, which is generi- 
cally impossible as shown above. 



2. The Cooperon propagator 

Finally, we return to the Cooperon propagator, which 
in our language is given by the correlation function 

Cn,n2,n,n, = j dy { (? (AQ)„^„^ (y) 

x?(A0)„3n.W)},j^ . (3-22) 
The decomposition analogous to Eq. (3.17) reads in this 



(7 = A A r7(dc) 

-1-712 ,n3-F"4 ^mri2,n3Ti4 ' (3.23) 

From Eq. (3.6b) we see that the disconnected piece can 
be expressed in terms of Green functions as 

= 1^ E {^k.p(^'^") ^-k,-p(i'^m)}dis • (3-24a) 



k,p 



Here Q is the electron Green function before disorder av- 
eraging. It is related to G, Eq. (1.1), by {5k,p(iw„)}dis = 
(5k, p G(k, iujn). Analogously, the diffuson propagator can 
be written 

=1^E{^k,p(ia;„)ap,k(ia;„)},i, . (3.24b) 



k,p 
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In the presence of time reversal invariance one has the 
identity 

Gk,p{iuJn) = G-p-kiiuJn) ■ (3.25) 
Substituting this into Eq. (3.24a) leads to 

Cit^ = -^it^ , (3.26) 

in time reversal invariant systems. This proves that the 
Cooperon is massless as well. 

IV. CONCLUSION 

In summary, we have considered a model of electrons 

in the presence of quenched disorder. For this model, 
we have proved that the diffuson propagator in d > 2 at 
T = remains massless in the presence of an electron- 
electron interaction. This means that the dephasing rate 
due interactions vanishes at T = 0. In systems that are 
time-reversal invariant, the Cooperon propagator is also 
massless. These results prove that, to the extent that 
time reversal invariance holds for the systems studied in 
Ref. 10, the calculations and argmncnts given in Rcfs. 13, 
14 cannot be used to explain the experimental observa- 
tions. 
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